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GeV emission from Gamma Ray Bursts: a radiative fireball? 
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ABSTRACT 

We study the emission observed at energies > 100 MeV of 1 1 Gamma Ray Bursts (GRBs) de- 
tected by the Fermi/Large Area Telescope (LAT) until October 2009. The GeV emission has 
three main properties: (i) its duration is often longer than the duration of the softer emission 
detected by the Gamma Burst Monitor (GBM) onboard Fermi [this confirms earlier results 
from the Energetic Gamma-Ray Experiment Telescope (EGRET)]; (ii) its spectrum is consis- 
tent with F v oc v~ l and does not show strong spectral evolution; (iii) for the brightest bursts, 
the flux detected by the LAT decays as a power law with a typical slope: t -1-5 . We argue that 
the observed >0.1 GeV flux can be interpreted as afterglow emission shortly following the 
start of the prompt phase emission as seen at smaller frequencies. The decay slope is what 
expected if the fireball emission is produced in the radiative regime, i.e. all dissipated energy 
is radiated away. We also argue that the detectability in the GeV energy range depends on the 
bulkLorentz factor T of the bursts, being strongly favoured in the case of large T. This implies 
that the fraction of bursts detected at high energies corresponds to the fraction of bursts hav- 
ing the largest T. The radiative interpretation can help to explain why the observed X-ray and 
optical afterglow energetics are much smaller than the energetics emitted during the prompt 
phase, despite the fact that the collision with the external medium should be more efficient 
than internal shocks in producing the radiation we see. 
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1 INTRODUCTION 

The Fermi Gamma Ray Space Telescope (Fermi) has onboard two 
instruments: the Large Area Telescope (LAT), sensitive in the 100 
MeV - 100 GeV energy range (and even beyond 100 GeV, for 
very bright sources, Atwood et al. 2009), and the Gamma Bursts 
Monitor (GBM), especially designed for the detection of Gamma 
Ray Bursts (GRBs), sensitive in the 8 keV - 40 MeV energy range 
(Meegan et al. 2009). The LAT revealed 12 GRBs above 100 MeV 
confirming that GRBs can be sources of very high energy photons 
and that the fraction of GRBs that can be detected at these energies 
is roughly 10 per cent of those detected by the GBM at lower ener- 
gies. It was the EGRET instrument, onboard the Compton Gamma 
Ray Observatory (CGRO) the first to detect GRBs above 100 MeV 
(Fishman & Meegan 1995; Kaneko et al. 2008), but it is the much 
better sensitivity (and reduced dead time) of the LAT to allow us 
for the first time to try to understand the origin of this emission 
and to answer the question: does it belong to the prompt phase or 
is it afterglow emission produced by the fireball colliding with the 
circum-burst medium? Or has it still another origin? 

One of the puzzling features of the high energy emission as re- 
vealed by EGRET was that it was long lasting, yet it started during 
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the prompt phase as seen by the Burst Alert and Transient Exper- 
iment (BATSE) onboard CGRO sensitive in the 30 keV - 1 MeV 
energy band. For instance, GRB 940217 emitted > 100 MeV pho- 
tons up to 1.5 hours after the prompt phase ended in the BATSE 
detector. A photon of 18 GeV was received ~5000 s after the trig- 
ger (Hurley et al. 1994), and this was the highest photon energy of 
a GRB until the Ferrai-LAT era. On the other hand, about a third 
of the high energy photons were received within 120 s, before the 
end of the prompt phase as detected by BATSE. 

Up to now, there have been three LAT-detected GRBs already 
discussed in the literature. In GRB 0809 16C (Abdo et al. 2009a), 
there is evidence that the spectrum from 8 keV to 10 GeV can be 
described by the same Band function (i.e. two smoothly connected 
power laws), suggesting that the LAT flux has the same origin of the 
low energy flux. On the other hand, the flux level of the LAT emis- 
sion, its spectrum and its long lasting nature match the expectations 
from a forward shock, leading Kumar & Barniol-Duran (2009) to 
prefer the "standard afterglow" interpretation (see also Razzaque, 
Dermer & Finke 2009 for an hadronic model; Zhang & Peer 2009 
for a magnetically dominated fireball model and Zou et al. 2009 for 
a synchrotron self-Compton origin). 

In the short bursts GRB 090510 the spectrum in the LAT en- 
ergy range is not the extrapolation of the flux from lower energies, 
but is harder, leading Abdo et al. (2009b) to propose a synchrotron 
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Table 1. The 12 bursts detected by the Fermi-LAT instrument above 100 MeV, until October 03 2009. Besides their redshifts (when measured) and duration, 
we give the parameters of the time integrated GBM spectrum collected from the literature and the corresponding reference. Fluences S are in [erg cm -2 ], 
peak energies -E pea]l in keV. In column 10 we report the fluence S in the [8 keV-10 MeV] energy range calculated from the spectral parameters of the GBM. 
Column 11 reports the fluence in the [0.1-100 GeV] energy range obtained from the analysis of the LAT spectra performed in this paper (whose results are 
given in Tab. 2). We adoped a BAND model for the GBM, and a simple power law of photon slope T for the LAT. When /3 is not indicated, the adopted fitting 
model is a cut off power-law of photon slope a. a : Sgbm in the [8 keV-30 MeV] energy range b : Sgbm m the [50 -300 keV] energy range; c : Sgbm m me 
[50 keV-40 MeV] energy range; d : Sgbm in the [50 keV-10 MeV] energy range. The number quoted in the "Ref." column refer to GCN circulars as follows: 
8141: van der Horst & Connaughton 2008; 8278: van der Horst & Goldstein 2008; 8407: Omodei 2008; 8682: Chaplin, van der Horst & Preece 2008; 8902: 
von Kienlin 2009; 9021: Ohno M. et al. 2009; 9057: Rau, Connaughton & Briggs 2009; 9336: Guiriec, Connaughton & Briggs 2009; 9579: von Kienlin 2009; 
9866: Bissaldi & Connaughton 2009; 9933: Bissaldi 2009; 9983: Rau 2009. 
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self-Compton interpretation for its origin. Instead we (Ghirlanda, 
Ghisellini & Nava 2009) proposed that the LAT flux is afterglow 
synchrotron emission, on the basis of its time profile and spectrum 
(see also Gao et al. 2009; De Pasquale et al. 2009). 

Finally, the LAT flux of GRB 090902B decays as t~ 1 ' 6 (Abdo 
et al. 2009c), it lasts longer than the flux detected by the GBM, and 
its spectrum is harder than the extrapolation from lower frequen- 
cies, making it a good candidate for an afterglow interpretation, de- 
spite the arguments against put forward by Abdo et al. (2009c), that 
we will discuss in this paper. Moreover, in GRB 090902B there is 
evidence of a soft excess (observed in the GBM spectrum below 50 
keV) which is spectrally consistent with the extrapolation at these 
energies of the LAT spectrum. 

As the few examples above demonstrate, there is no consen- 
sus yet on the nature of the high energy emission of GRBs. Since 
only three of the nearly dozen bursts detected by the LAT have al- 
ready been discussed in the literature, we present here a study of 
the entire sample of bursts detected at high energies by the LAT. 
We will construct the light curves of the high energy flux and the 
spectral shape in the 0.1-100 GeV energy range, to find if there are 
properties that are common among different bursts that can help to 
understand their nature. 

Indeed, we believe that a consistent scenario emerges: the LAT 
spectra are often inconsistent with the extrapolation of the GBM 
spectra (except two cases) and the light curves can be often de- 
scribed by a power law decay in time, i.e. Flat oc t~~ a , with a 
slope often close to a = 1.5. In the brightest cases also the ris- 
ing part is visible, and is consistent with Flat oc t 2 . These are, in 
our opinion, strong indications of the afterglow nature of the LAT 
emission. Furthermore, we suggest that GRBs with a flux decay- 
ing as Flat oc t~ , and with a spectral slope around unity [i.e. 
F{v) oc v~ ], are emitting in the radiative regime of a forward 
shock. We will also point out the role that the electron-positron 
pair production process has in establishing the radiative regime. Fi- 
nally, we will discuss the consequences of our findings. 



We adopt a cosmology with h = S1a = 0.7 and Qyi — 0.3 
and the convention Q = 10 X Q X , using cgs units. 



2 SAMPLE AND DATA ANALYSIS 

We considered all the 12 bursts detected in the Fermi-LAT until 
the 3rd of October 2009. These are reported in Tab. 1 with their 
redshifts (Col. 2) and the spectral parameters and fluences (Col. 4, 
5, 6, 7) as reported in the literature, obtained from the analysis of 
the GBM spectrum. Since the GBM fluences reported in the lit- 
erature refer to different energy ranges, we convert all the GBM 
fluences to the common 8 keV - 10 MeV energy range (Col. 10). 
In addition (last column) we report the fluences in the 0. 1-100 GeV 
energy range of the LAT obtained from the spectral analysis of the 
LAT data (spectral parameters are given in Tab. 2). For those GRBs 
with measured redshifts we computed the isotropic equivalent en- 
ergy F 7 ,i so by integrating the GBM spectral model in the in the 1 
keV - 10 MeV rest frame energy range. 

Among the considered bursts there are three cases which 
have been published in recent papers: GRB 0809 16C (Abdo et al. 
2009a), GRB 090510 (Abdo et al. 2009b; Ghirlanda et al. 2009) 
and GRB 090902B (Abdo et al. 2009c). All the others are unpub- 
lished. We did not consider GRB 081215 which, lying at a large 
angle (86°) with respect to the LAT boresight (Preece et al. 2008), 
required a non standard analysis of the LAT data (McEnery et al. 
2008). GRB 081024B and GRB 090510 are of the short class. 
Seven bursts have measured redshifts, for all the others we assume 
a typical redshift of 2 and 1 for the long and short class. 

We have analysed the Fermi-LAT dauQ with the Fermi 
ScienceTools (v9rl5p2) released on Aug. 8th 2009. LAT 
count light curves (extracted with the gtbin tool) were rebinned 
in time with a variable bin size, different for each burst. 

1 http://fermi . gsfc . nasa. gov/ssc/data/ 
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Figure 1. Light curves of the 11 GRBs detected by LAT plus GRB 940217, as detected by EGRET (bottom right panel). The hatched region represents the 
duration (T90) of the emission detected by the GBM in the 8 keV-40 MeV energy range (for GRB 940217 it refers to the emission detected by BATSE). Times 
are in the observer frame for all bursts and arrows represent 2<r upper limits. 



We analysed the spectrum of the emission detected by the 
LAT. For the brightest part of the burst we applied the standard 
procedure (i.e. extracted the spectra and created the relative re- 
sponse files with the gtbin and gtrspgen tools, respectively). 
We considered the spectrum over a time interval covering entire 
light curve, and if the burst was particularly bright we also extracted 
the spectrum over a time interval coincident with the duration of the 
emission as observed by the GBM. To verify if and at what extent 
the LAT spectrum could vary with time, we extracted a series of 
consecutive spectra for each burst. As in most bursts we did not 



find evidence for substantial spectral evolution of the LAT compo- 
nent, we used the average spectrum to convert the count rate into 
physical units. 



3 RESULTS 

Light curves - Fig. Q] shows the light curves obtained from the se- 
lection of the LAT events with energies > 0.1 GeV. In each plot 
we also show the time interval (hatched region) corresponding to 
the duration of the GBM light curve (T90 in Tab. 1). In 9/1 1 events 
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there is a peak in the LAT light curve and the latter has a duration 
much longer than the duration of the GBM light curve (shown by 
the hatched region in Fig.QJ. After the peak, the light curves of dif- 
ferent GRBs show a similar temporal decay. In a few cases (see also 
Ghirlanda et al. 2009) a rising of the light curve as t 2 is seen before 
the peak. The three faintest GRBs (GRB 090323, GRB 090328 and 
GRB 090626) have light-curves that appear much flatter than the 
other ones (please note the different scale of their y-axis) and we 
cannot exclude that the background, in this cases, plays some role. 
The bottom right panel shows the light-curve of GRB 940217 as 
detected by EGRET (Hurley et al. 1994), selecting photons above 
100 MeV. As can be seen, also this burst show a similar decaying 
light curve. 

Spectral evolution - In Tab. 2 we report the results of the LAT 
spectral analysis. For each burst the first line refers to the spectrum 
used to convert the count rate into physical units while the follow- 
ing lines give the spectral index for each time resolved spectrum. 
We report in Tab. 2 also the flux integrated between 100 MeV and 
100 GeV. By comparing the time resolved spectral results of in- 
dividual bursts we see that there is no evidence of strong spectral 
evolution of the LAT spectral index during the burst. On average, 
all the spectral index are distributed between 1.5 and 2.2. 

Spectral slopes in the LAT vs GBM - In Fig. [2] we compare the 
spectral index of the LAT emission with the spectral index of the 
average GBM spectrum (whose spectral parameters are reported in 
Tab. 1). The low energy spectral index a (circles in Fig. [2] red in 
the electronic version) of the Band model (or of the cutoff power- 
law model for GRB 081024B and GRB 090323) is systematically 
harder than the spectral index of the LAT component. The high en- 
ergy spectral index /? of the Band model (open squares in Fig. [2] 
blue in the electronic version) appears softer than the LAT spec- 
trum. An extreme case is GRB 090902B which clearly shows that 
the LAT component is spectrally different from the tail of the Band 
function. Indeed, in this burst there is also evidence of a soft spec- 
tral excess detected in the GBM below 50 keV (Abdo et al. 2009c; 
De Palma et al. 2009). We also note that in only two bursts, GRB 
080916C (Abdo et al. 2008) and GRB 090926 the high energy spec- 
trum of the Band model is consistent with the spectral slope of the 
LAT data. 

LAT vs GBM fluences - Fig. [3] shows the fluence in the LAT en- 
ergy range 100 MeV - 100 GeV (using the fluxes listed in Tab. 2) 
as a function of the fluence in the GBM energy range 8 keV - 10 
MeV. The shaded regions correspond to 1, 2 and 3 a of the distri- 
bution of GBM fluences for the 121 GRBs detected so far by the 
GBM with measured prompt phase emission peak energy (Nava et 
al., in preparation) and that appeared in the Gamma Ray Bursts Co- 
ordinate Network (GCN) circularfl The dashed line marks equal- 
ity between the two fluences. We can see that all but the two short 
bursts (GRB 081024B and GRB 0905 10) have GBM fluences much 
brighter than average. If all GRBs with GBM fluences la brighter 
than average and in the LAT field of view (i.e. one half) were 
detected by the LAT, we should have a fraction of LAT-detected 
GRBs of ~ 16%, that is not far from the actual fraction (see also 
Guetta & Pian 2009). One can compare Fig.|3]with Fig. 4 of Le & 
Dermer (2009), showing the pve-Fermi bursts detected by EGRET 
and BATSE. Apart from GRB 930131, showing an EGRET fluence 
comparable to the BATSE one, all the other pre-Fermi GRBs seem 

2 http://gcn.gsfc.nasa.gov/ 
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Figure 2. Spectral index (a and j3, red circles and blue squares respec- 
tively) of the GBM time integrated spectra (reported in Tab. 1) versus the 
spectral index obtained from the analysis of the LAT data presented in this 
paper (Tab. 2). The dashed line represents equality. The two lower limits are 
GRB 081024B and GRB 090323 whose time integrated GBM spectrum is 
best fit by a cutoff power-law model. For illustrative purposes we assumed 
fi > 4 for these two bursts. This plot shows that the LAT spectrum is softer 
than the low energy spectral index of the Band model fitting the GBM spec- 
trum (red circles) and it is harder than the high energy spectrum of the Band 
model fitting the GBM spectrum (blue squares). 

to be characterized by a fainter GeV fluence relative to their flu- 
ence at smaller energies, but the sample is too small to draw any 
conclusion. 

Time decay of the LAT flux - We converted the count rate of 
Fig. [3 into luminosity. For the bursts without measured redshifts 
we assumed a typical redshift of 2 for long events, while for GRB 
081024B we used a redshift of 1. We show the light curves of 8 
GRBs with good quality data in the top panel of Fig. [4] where the 
times are in the source rest frame. The grey shaded stripe has a 
slope of t~ w / 7 , and it is shown for comparison. We can see that 
the light-curves show a power-law behaviour, and that the decay 
slope is often steeper than unity. Initially, some bursts show a ris- 
ing phase and therefore it is possible to define the peak time of 
their high energy emission. As seen below, if the peak time marks 
the onset of the afterglow emission it can be used to estimate the 
bulk Lorentz factor T. 

Common decay for the brightest LAT bursts - The bottom panel 
of Fig. |4]shows the light curves of the 4 brightest GRBs with red- 
shift, once the 0. 1-100 GeV luminosity is divided by the energetics 
Ef,iso of the flux detected by the GBM. The shaded stripe has a 
slope t~ 10 / 7 , and it is shown for comparison. These four GRBs 
show a common behaviour, being all consistent, within the errors, 
with the same decay, both in slope and in normalisation. Note that 
GRB 090510, a short bursts, behaves similarly to the other 3 bursts, 
that belong to the long class, but its light-curve begins much earlier. 
If we divide the light-curves by the average luminosities as derived 
by the GBM [instead of the energetics; i.e. we multiply by the time 
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Table 2. LAT spectral results. We give the time interval (to — <i) for the ac- 
cumulation of each spectrum, the photon index, the flux integrated between 
100 MeV and 100 GeV, and the C-statistic for degrees of freedom. Errors 
are given at the 90% confidence level. 
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290/(1+2:)] the resulting light-curves of the 4 GRBs spread within 
a larger region. 

To conclude, we find that i) the LAT fluxes decay as a power- 
laws; ii) the spectral shape at high energies is not strongly evolving; 
iii) the LAT spectrum has a slope intermediate between the low and 
high energy slope (i.e. a and f3) of the Band function used to fit the 
GBM data; iv) the brightest 4 GRBs show a common t — 1,5 decay 
and even the same normalisation, once their LAT luminosities are 
divided by the GBM energetics. 

These characteristics are the same as observed/predicted by 
the external shock scenario giving rise to the afterglow. We there- 
fore suggest that the high energy emission of the GRBs detected by 
the LAT has an afterglow origin. The fact that the high energy emis- 
sion overlaps in time with the prompt phase as seen in the GBM can 
be explained by invoking a relatively large value of the bulk Lorentz 
factor, corresponding to relatively small deceleration radii and on- 
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Figure 3. Fluence in the [0.1-100 GeV] LAT energy range as a function 
of the [8 keV - 10 MeV] GBM ones. Short GRBs are marked with filled 
squares, long GRBs with filled circles. GRBs with known redshifts are the 
ones with a LAT fluence larger than 2 X 10~ 5 erg cm -2 (red in the elec- 
tronic version). The shaded areas indicate the 1-2-3 cr values of the distri- 
bution of fluences of the 121 GRB with £ pcak (as of Oct. 2009) detected 
by the GBM. 

set times largely contracted by the Doppler effect. What is at odd 
with respect to the "standard afterglow" scenario is the relatively 
steep slope of the flux decay, even when the high energy spectrum 
indicates that we are observing this component close to its spectral 
peak. We offer a solution to this problem in the next section, where 
we will also argue that the likely emission process producing the 
high energy flux is synchrotron radiation. 



4 THE BOLOMETRIC AFTERGLOW LUMINOSITY 

In the early afterglow phases, the emission is likely to occur in the 
fast cooling regime, in which all the energy of the accelerated elec- 
trons is radiated away. In this case the bolometric afterglow lumi- 
nosity can be calculated in a simple way. Assume that the shock 
generated by the fireball has reached a radius R, and that it moves 
within a region characterised by a uniform number density n (this 
case can be easily generalised to different density radial profiles). 
The (comoving) emitting volume is V' = A-kR 2 AR', since we are 
assuming that the fireball is a spherical shell. The radiative cooling 
rate of the electrons is measured by 7 where 7m c c 2 is the electron 
energy, and the emitting particles are distributed in energy accord- 
ing to N(y). Note that the time derivative, the electron energies and 
their energy distribution are all measured in the comoving frame. 
In this case the bolometric luminosity is: 



V'N(-y)-yd>y 
4nR 2 F 2 m c c 2 J N{-y)-yAR'{~/)d>y 



4nR 2 r 2 m c c 3 I ^(7)7^7 



(1) 
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Figure 4. Top panel: Light curves of the 8 brightest bursts GRBs detected 
by LAT. The luminosities are integrated in the 0.1-100 GeV energy range 
at the source rest frame. For GRBs without measured redshifts we assumed 
z = 1 for short and 2 = 2 for long events. The time is in the rest frame of 
the sources. Upper limits are at la level. The grey stripe indicates a slope 
j— 10/7 ^ B tt om panel: Light curves of the 4 brightest GRBs with redshift, 
normalised to the total energetics of the GBM energetics. The luminosities 
are integrated in the 100 MeV-100 GeV energy range at the source rest 
frame. For GRBs without measured redshifts it is assumes 1 for short and 
2 for long events. The time is in the rest frame of the sources. Upper limits 
are at la level. The yellow stripe indicates a slope t~ 10 / 7 . 



We have used the fact that the distance AR' can be approxi- 
mated by the cooling length as measured in the comoving frame: 
AR'("f) — ct' cool — C7/7. Therefore AR' is energy dependent, it 
is smaller for high energy particles, that spend most of their energy 
faster. Eq. Q] is remarkably independent of the specific radiation 
process. The integral in Eq.[T]must correspond to the fraction e c of 



the available energy density as measured in the comoving frame, 
i.e. 



(2) 



m c c / N^jdj — e c nT m p c 

Therefore Eq.Q]becomes: 

Liso = 4nR 2 F 4 m p c 3 e e n 

= 16ira 2 t 2 T 8 m p c 5 e e n 



(3) 



where we have assumed that the size R is measured by the observed 
time as R = 2actT 2 . The factor a is equal to 1 if the fireball moves 
at a constant speed, and becomes greater than 1 when it decelerates 
(see e.g. Sari 1997). Eq.|3]is valid as long as the afterglow is in the 
fast cooling regime, irrespective of the radiative or adiabatic nature 
of the process, that changes only the relation between the observed 
time t and the bulk Lorentz factor F at that time. In fact, when the 
forward shock is coasting (i.e. before being notably decelerated) 
we have Li so oc t 2 in both cases. When the shock starts to deceler- 
ate, the observed luminosity decreases according to the appropriate 
r(t) function, that is different for the adiabatic and radiative cases. 

Adiabatic case - We adopt the following relation between the ob- 
served time and V: 



3E k 



32-Ka 3 nm p c 5 t 3 



(4) 



where £?k,f is the kinetic energy of the fireball after the prompt 
phase. The same equation can be used to define the deceleration 
time idec, once we set a — 1 and substitute To to F. If r\ is the effi- 
ciency of conversion of the initial kinetic energy i?k,o into radiation 
of the prompt phase, we have 



-Ek.f = E\r 



l-T) 



(5) 



When the fireball is still in its coasting phase the observed luminos- 
ity increases as t due to the increased visible area. After t^ ec the 
observed luminosity decreases as t~ , as can be seen inserting Eq. 
HintoEq.[U 



-^iso,i 
Liso,; 

^peak,; 



l6iYt 2 rlrn p c 5 e B n; t < i pcak 

3 tc-Ek.f 
2a 



t 

3-Ek,f 



327rarQ?impC 5 



3£ k ,f 
32irFgnm p c 5 



1/3 



1/3 



t > t p 

^dec 



(6) 



To find fpeak we equated the two expressions for Li s 



Radiative case - In this case an important fraction of the dissipated 
energy is radiated away. This implies that the emitters, i.e. the elec- 
trons, receive a large fraction of the available energy (directly or 
through the interactions with protons, and/or through reconnection 
of the magnetic field) and radiate it efficiently. In this case the en- 
ergy of the fireball decreases, changing the T(t) function. This has 
been studied by Blandford & McKee (1976); and the solution is 
(Katz & Pian 1997; Vietri 1997; Sari, Piran & Narayan 1998): 



r 

x 



(r + i)(x + i) 2 + (r -i) 
(r + i)(x + i) 2 -(r -i) 



m 



4-kF m p ncR 3 
3-Ek.f 



(7) 
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where Mf is the mass of the fireball and m is the swept interstellar 
mass. When the fireball is decelerating, but still relativistic, X C 1 
and Eq.|7] simplifies to: 



3-Ek.f 



1/7 



,3/7 
T> dec _/. 



-3/7 



(8) 



Inserting this into Eq.[3]we obtain: 



^peak,r — 



167rt 2 ro»7ipC 5 e e n; 

3tc-Ek,f 3/7 ,-10/7 
2a 10/7 Mec 

^dcc 



oak 



5/12 



(9) 



The peak time of the bolometric afterglow emission (estimated 
equating the two limiting forms of Li so ) precedes the decelera- 
tion time by a small factor. Integrating dR = 2cT 2 dt assuming 
r oc t~ 3 ^ 8 (adiabatic) or V oc r.~ 3//7 (radiative) we have a = 4 or 
a = 7 for the adiabatic and radiative case, respectively. Therefore 
tpcak = 0.63tdcc (adiabatic) and t pcak = 0.44£ dcc (radiative). 

After the peak time, radiative afterglows decrease faster than 
adiabatic ones, as the fireball energy is no longer constant, but de- 
creases. As noted by Sari, Piran & Narayan (1998), partially ra- 
diative fireballs would have scalings intermediate between the pure 
adiabatic and pure radiative limits. Even if, initially, a fireball is 
purely radiative, after some time it must become adiabatic, as a con- 
sequence of incomplete cooling of the accelerated electrons. If the 
electrons are accelerated above some minimum energy 7 m m e c 2 , 
this will occur when this electrons cannot cool in a dynamical time, 
so when j m = j c , where 7 c m e c 2 is the energy of those electrons 
cooling in t' ~ Tt ~ R/(acF). 

When observing the flux in a particular frequency range Av, 
we are never observing the bolometric flux, so in general the time 
decays are different from t -1 (adiabatic) or t~ 10 ^ 7 (radiative). If 
the emitted spectrum (in a vF v plot) has a peak at ^ pca k, and i/ m 
decreases in time, then the time decay would be flatter for v < 
^pcak, and steeper for v > ^ pca k- However, if the observed flux has 
a spectral index close to unity (i.e. v ~ f P eak)> then the observed 
flux becomes a good proxy for the bolometric one, with the same 
time decay slope. 

For a uniform circum-burst medium, the relation between the 
decay slope a and the spectral index /3 for a flux density F(y, t) oc 



7 



is (Sari, Piran & Narayan 1998): 

(6/3 - 1) 



(10) 



returning a — 10/7 = 1.43 when (3 = 1 and a = 1.77 for 
P — 1.2. This derivation assumes that the number of accelerated 
electrons is always a fixed fraction of the protons present in the 
circum-burst medium. 



5 PAIR-ENRICHED INTERSTELLAR MEDIUM 

When the prompt phase emission spectrum extends above 
Spe a k(l + z) ~ m e c 2 we can convert a fraction of the high energy 
photons into electrons-positron pairs. This case has been studied in 
detail by Thompson & Madau (2000), Meszaros, Ramirez-Ruiz & 
Rees (2001) and especially by Beloborodov (2002). 

The basic idea is that although the scattering depth of the cir- 
cumburst medium can be much smaller than unity, it can neverthe- 
less scatter a fraction of the prompt phase photons along non ra- 
dial directions. These scattered photons can then interact with the 



arriving high energy prompt phase photons producing pairs. The 
process is not controlled by the probability of the interaction be- 
tween the scattered and the primary prompt phase photons: this is 
almost unity (up to very large distances), due to the huge amount 
of the prompt phase photons. The process is controlled by how 
many photons are scattered. The full description of this scenario is 
rather complex, and we refer to Beloborodov (2002) for the com- 
plete treatment. We focus here on a few estimates, to give the idea 
of the importance of the process. The basic quantity of interest is 
the number of scatterings done by a single electron located at a dis- 
tance R from the emission site of the prompt phase emission. Using 
the Thomson cross section for simplicity, and setting hv = xm e c 2 , 
this number is 



(XT 



{x}m e c 2 A7YR 2 ct huT; 



||t0 -E 7 ,iso,54 (n) 



Almost all these photons will be converted into pairs immediately 
after they have been scattered. This implies that the circumburst 
medium will be greatly enriched by pairs before the arrival of the 
forward shock. This can occur even if the total number of the in- 
tercepted photons is a tiny fraction of the total. For instance, if 
the interstellar medium is homogeneous with density n, the to- 
tal number of scattered photons within 10 17 cm is only a fraction 
tt = 6.65 x 10~ 8 n of the total number of photons of the prompt 
phase. But this is enough to greatly pair-enrich the circumburst 
medium. Furthermore, the scattering and the pair production pro- 
cesses pre-accelerate the interstellar medium. If there is one pro- 
ton per primary electron, and if the energy deposited by the sin- 
gle scattering with subsequent pair production is roughly equal to 
m c c 2 , this process will be important below a certain distance, be- 
low which there occur more than 1000 scatterings for primary elec- 
tron (i.e. in this case the proton associated with the primary electron 
will start to move with V ~ 2 in the radial direction). As a feed- 
back, if the medium starts to move then the typical energy of the 
scattered photons will start to decrease, quenching off the pair pro- 
duction process (i.e. the scattered photons have too small energies 
to interact with photons around a few MeV). On the other hand, the 
produced pairs, if are re-isotropized in a short time, can also scat- 
ter the incoming prompt phase radiation, enhancing the process and 
making it exponential. 

Therefore Eq.QTJis only a simple but rough estimate of a much 
more complex scenario. We can nevertheless draw some important 
conclusions: 

• Pairs are important if the prompt phase emission extends 
above threshold. 

• At a negligible expense (i.e. the fraction of absorbed prompt 
phase emission is negligible) the environment is largely enriched 
by pairs. 

• The distance for which the number of produced pairs equals 
the number of primary electrons is sufficiently large and affects the 
properties of the forward shock up to some relevant observed time. 
For instance, the "closure" relation given by Eq.[l0]is modified as 
long as the number of pairs per proton is larger than unity, because 
in this case the energy 7 m oc Yn/n+ oc FR 2 . Here n+ is the pair 
density. Introducing this extra R 2 dependence we find 

a = |(4/3 + 1) 



10/7 = 1.43 when /3 = 1 and a 



(12) 
1.66 for 



returning a 
P = 1.2. 

• Although the details of the shock acceleration process are con- 
troversial, it is reasonable to assume that the ratio of the energy 
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given to leptons and protons will increase, if we have many leptons 
per proton. This is then one way to have a radiative fireball. 

We therefore propose that bursts whose prompt phase emission ex- 
tends above m c c 2 should be characterised by an early radiative 
(then powerful) afterglow. 



5.1 Additional processes 

We consider here other processes that can be relevant for the for- 
mation of the high energy afterglow: 

• When fbursts > £dcc the region of the forward shock where 
leptons are accelerated is illuminated by the flux of the prompt 
phase emission (of luminosity L 7l i so ). This component lasts as long 
as the forward shock is illuminated by the prompt phase (see Be- 
loborodov 2005a). The corresponding energy density, as measured 
in the comoving frame of the forward shock is 



(13) 



A-kR?cT 2 

where the subscript "ext" stands for "external" to the afterglow 
emitting region. This has to be compared with the local magnetic 
energy density 



t/ B = <LB,nm v cT 2 



(14) 



Therefore the ratio between the synchrotron and the "external 
Compton" (i.e. the luminosity produced by scattering U^) lumi- 
nosities is (see also Beloborodov 2005a): 



Lec 
L s 



' cxt 



47ri? 2 r 4 e B ? 



= 0.18 



fL 



7,iso,53 



Rl 7 Tje B ,- 



■(15) 



The factor / < 1 accounts for the suppression of the power emitted 
in the direction of the observer due to the anisotropic pattern of the 
incoming photons in the frame of the fireball. An order of magni- 
tude estimate of its value can be gained through a simple example. 
In the frame of the fireball, assume that all the seed photons for the 
scattering are coming radially. Electrons travelling at 8' = 180° 
from the photons lose energy at a rate oc j 2 (l — [3 cos 9') ~ 47 2 . 
Electrons moving at 90° lose energy at a rate oc j 2 . This is the 
emission that the observer (on the Earth) will preferentially see. 
Therefore the factor / is less than, but of order of, unity. This ex- 
ternal Compton component would start to be important at frequen- 
cies above u E c ~ Tm^peak ~ 7m,3^ P cak,M e v TeV. Below u E c we 
should have F(v) oc ^~ 1//2 . 

• The high energy emission can also be produced by the syn- 
chrotron self-Compton (SSC) process (See e.g. Corsi et al. 2009; 
Fan et al. 2008), particularly important when i) e c > 6b ; ii) we are 
in the fast cooling regime and iii) we are in the Thomson limit (i.e. 
the scattering can be described by the Thomson cross section). 

Condition i) and ii) are always fulfilled in radiative fireballs, 
while condition iii) may be violated. The limit for the Thomson 
regime can be derived considering the dimensionless frequency 
x' m — hvm/(rm e c 2 ) (as measured in the comoving frame) and 
the electron energy 7 m . If x' ln 'y m > 1 The entire process occurs in 
the Klein Nishina regime if x m 7m > 1, i.e. when: 



■p3 3 1/2 1/2 



> 1.77 x 10 1 



(16) 



For moderate pair production (i.e. n+/n< 100) and for still large 
r the early SSC process is then in the Klein Nishina regime, and 
is therefore inefficient. Furthermore, the SSC spectrum starts to be 
important at z^ssc given by 



7m 



n .. 1r ,22 4 1/2 1/2 

7 X 10 e„e n n 1 



Hz (17) 



It is a strong function of n/n+: for less than 100 pairs per proton 
(and still a large F) the SSC spectrum starts at frequencies above 
the LAT range (with a flux reduced by Klein — Nishina effects). 
The mid panel of Fig. [5] shows i^ssc as a function of time for one 
particular case. 

We conclude that the most likely radiation process originating 
the LAT emission is synchrotron. 

To illustrate the above considerations and to give an example 
of the predicted high energy flux in radiative fireballs, we have cal- 
culated the bolometric flux emitted in one specific case, assuming 
that the prompt phase energetics -E 7 ,i so = 10 53 ergs, r\ — 0.2, 
z = 1, To = 1000, n = 1 cm -3 , p — 2. Furthermore, we 
assumed a duration of 1 s and e c = 0.9, eb = 0.1. The re- 
sulting bolometric luminosity (normalised to _E 7; i so ) is shown in 
the top panel of Fig. [5] together with its corresponding energet- 
ics [i?bol(i) = J Lboi{t')dt']. We also indicate the t _1 and 
the t~ 10 / 7 time behaviour (dashed black lines). The mid panel 
shows the time profile of 3 characteristic frequencies: the injected 
frequency v m , the cooling frequency v c and the SSC frequency 
fssc = 7^111 (see also Beloborodov 2005b for the case of pair 
enriched circum-burst material, but with an adiabatic fireball). 

The 2 upper shaded areas correspond to the frequency ranges 
covered by the LAT and GBM, while the lower one indicates the 
optical frequency range. The bottom panel shows the time profile of 
the minimum Lorentz factor of the injected electron 7 m , the cool- 
ing Lorentz factor y c , the bulk Lorentz T, together with the time 
profile of the magnetic field B and the number of pairs per pro- 
ton n+ /n, calculated according to Eq.[TT] This quantity is crucial 
to calculate 7 m , since the same available energy must be divided 
by the total number of leptons, including the pairs. Since their 
amount changes with R (and correspondingly with the observed 
time), the time profile of 7 m is greatly modified by the presence 
of pairs. As a consequence, both Vm and fssc are largely affected, 
their values being much lower than in the absence of pairs. A note 
a caution: although the presence of pairs may be crucial to bring 
the process to the radiative regime, the exact amount of pairs is 
difficult to calculate, being partially dependent on the exact shape 
and time evolution of the spectrum of the prompt phase emission 
above threshold, the presence or not of a magnetic field embedded 
in the circum-burst medium, a possible clumping of this medium, 
and so on. Ours are bound to be only rough estimates. Bearing 
the above caveat in mind, we find that the synchrotron emission, 
at the peak time, should have a flux F(y) oc z/ -0 ' 5 between v c 
and Vm and F{v) oc v~ v ^ 2 (equal to u~ x in this example) up to 
^max — Vm (7max/7m) 2 - Therefore /max /7m>10 3 ensures that 
the synchrotron emission extend up to the GeV range. 

Note that the v~ p ^ 2 part of the spectrum may start in the GBM 
energy range, depending on the exact amount of pairs. There is 
then the possibility that the afterglow emission "contaminates" the 
prompt phase emission seen by the GBM. In some cases, this "con- 
tamination" can appear as an excess at both extremes of the GRB 
energy range, especially if pairs are very important, decreasing 7 m 
(as in the case of GRB 090902B, Abdo et al. 2009c). Also the op- 
posite (i.e. the prompt phase "contaminates" the afterglow seen in 
the LAT) can occur, especially when the high energy Band index f3 
is not too soft. In this latter case most of the prompt phase photons 
contributing to the LAT flux should be at low energies. 
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Figure 6. The long burst GRB 080916C. Parameters are listed in Tab.[5] 



Figure 5. Top panel: time profiles of the bolometric luminosities and the 
corresponding cumulative energetics, in units of the initial kinetic energy of 
the fireball. For this particular example, we have assumed a radiative fireball 
with z = 1, £ 7 , iso = fO 53 erg, r] = 0.2, T 90 = Is, T = 10 3 , e c = 0.9 
and p = 2. The circumburst medium is homogeneous with density n = 1 
cm -3 . The dashed lines corresponds to t^ 1 and t~ 10 / 7 j i e the adiabatic 
and radiative cases. Pair production is accounted for in a approximated way, 
assuming that all scattered photons are transformed into pairs, but assuming 
that there are at most m p /m c pairs per primary electron. Mid panel: the 
time profiles of the frequencies i/ m , u c and ^SSC = 7m* / m- The hatched 
areas mark the energy ranges of the LAT instrument [0.1-100 GeV], the 
GBM instrument [8-1000 keV] and the optical range (corresponding to the 
U and R filters). Bottom panel: the time profiles of the injected energy *y m 
and the cooling energy 7 C . We also show the profile of T, of the magnetic 
field B (assuming eb = 0.1), and the number of pairs per primary electron 
n+/n. Since n = 1, this also corresponds to the density of pairs. 



Table 3. Parameters for the radiative afterglow models 
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For simplicity, we have assumed that e c is constant, and not 
proportional to the amount of pairs per proton (since this number 
is uncertain). However, the radiative phase should end in any case 
when 7c becomes greater than 7 m since in this case most of the 
energy given to electrons cannot be radiated away in a dynamical 
time. 



6 APPLICATION TO SPECIFIC BURSTS 

We applied the radiative scenario to the 4 brightest (in the LAT) 
GRBs with redshift. They are the same illustrated in Fig. [4] namely 
GRB 080916C, GRB 090510, GRB 090902B and GRB 090926. 
In principle, the number of parameters used for the adopted model 
is limited (they are listed in Tab. [3), but we adopted a few rather 
drastic simplifications: 

• We consider the fireball, when colliding with the interstellar 
medium, as "thin". In other words, we assume that it can act as a 
piston having a total energy i?k,f • This is completely right for short 
GRBs, but not for long ones. According to Fig.[JJ the FermifLXT 
emission of several GRBs starts while the emission seen by the 
GBM has not ended. In this case the t 2 rising behaviour of the LAT 
light curve can be different (see Sari 1997). 

• When calculating the number of pairs produced by the 
circum-bursts medium, we neglect the amplification (exponential) 
effect of the produced pairs that can themselves scatter the in- 
coming radiation. The momentum deposited in the circum-bursts 
medium is also taken into account only by imposing that the maxi- 
mum number of pairs per proton is m p /m c , since a larger number 
corresponds to a mildly relativistic motion of the medium, and the 
quenching off of the pair-producing mechanism For simplicity, we 
use the Thomson cross section for scattering, and assume that most 
of the prompt phase photons are close to the threshold for pair pro- 
duction. 

• We assume that all electrons and positrons are accelerated. If, 
instead, only a fraction of them receive the entire available energy, 
then the typical Lorentz factors of the accelerated leptons is larger. 

• We use a fixed value of e c , even if the number of pairs popu- 
lating the circum-bursts medium decreases with R. Consequently, 
we use the radiative solution all throughout the shown evolution, 
with no transition to the adiabatic case. 
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Figure 7. The short burst GRB 090510 assuming T* = 0.6 s. Parameters 
are listed in Tab. [5] 



Bearing in mind these caveats, Figs.[6}{9]show the light-curves 
of the 4 GRBs interpreted on the basis of our radiative model, with 
the main parameters listed in Tab. [3] In the cases of GRB 0809 16C 
and GRB 090510 we have also added a constant flux to the light 
curve, to account for the presence of the background, flattening off 
the observed light-curves. In the case of GRB 0905 10, the fact that 
the flux above 200 s is due to the background has been confirmed 
by De Pasquale et al. (2009, see their Fig. 1). Also for 080916c the 
points above 1000 s are affected by background (see Abdo et al. 
2009a, and their Fig. 4). So, for these two bursts, the flattening of 
the light-curve at late times should not be due to the contribution of 
the SSC component entering in the LAT energy range (as predicted 
by Dermer, Chiang & Mitman 2000, and tantalisingly suggested 
by Fig. [6] and Fig. [7}, but only because we did not subtract the 
background. 

The solid lines shown in all top panels refer to the luminosity 
integrated in the 0. 1-100 GeV energy range, while the dashed thick 
lines are the bolometric fluxes (both normalised to the prompt phase 
energetics of each burst). For comparison we show also the lines 
corresponding to t" 1 and to t~ 10 ^ 7 . We can see that in all cases the 
radiative interpretation is in good agreement with what observed, 
and that in all cases the predicted v m is well below the 0.1 GeV 
value. This ensures that in the LAT we should see a spectral shape 
F(v) oc v~ p ^ 2 . The observed decay slope and the spectral index 
in the LAT energy range (see Tab. 1) are consistent with Eq.[l2] but 
the errors on ft — Tlat — 1 are too large to use this as a reliable 
test. 



7 DISCUSSION 

The found bulk Lorentz factors are in the range 630-900 for the 
long bursts, and 2000 for the short GRB 090510. We believe that 



Figure 8. The long burst GRB 090902. Parameters are listed in Tab.[5] 



these relatively large values are the key to understand why only a 
minority of bursts are detectable by the LAT. A large bulk Lorentz 
factor, in fact, means an early peak time of the afterglow (see Eq.[6] 
and Eq.|9j, and this in turn means a large flux. Faster fireballs have 
brighter afterglows. This is true for adiabatic as well as radiative 
fireballs. If the emission occurs in the radiative regime then the 
afterglow will be brighter still, since all the energy dissipated in the 
external shock is radiated away. 

If the circum-bursts medium is enriched by electron-positron 
pairs, we have a more favourable set up for a radiative process. If 
the acceleration mechanism divides its energy to all particles, then 
leptons should receive a total energy exceeding the one given to 
protons. But this may be only one of the means to have a radia- 
tive fireball. An alternative is to have a strong coupling between 
electrons and protons, with an efficient energy flow from protons to 
electrons. In any case, we can easily test if pairs are indeed im- 
portant by simply comparing the general properties of the early 
afterglow for bursts of different i5 pca k and high energy index /3, 
since only those bursts whose prompt phase photon energies ex- 
ceed m c c 2 should efficiently populate the circumburst medium by 
pairs. As an example, we may test if the high energy emission is 
present only in GRBs of high i5 pca k (in the rest frame) as it ap- 
pears to be the case until now, or if it occurs also for bursts with 
a small -E pca k- If this will occur, and if the flux will decay with a 
slower rate than t~ 10 / 7 , then we will have an indication of a fast 
fireball that emits adiabatically because of no pairs-enrichment of 
the circum-bursts medium. In other words, a possible test of the 
idea of having radiative afterglows because of pair enrichment is to 
find a different time decay for the high energy emission in classical 
GRBs whose prompt phase emission extends to high energies and 
X-ray flashes, characterised by relatively small values of S pea k. 

The radiative interpretation could ease the efficiency problem 
of the afterglow phase. This problem concerns the ratio of the ener- 
getics emitted during the prompt and afterglow phases, that is much 
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Figure 9. The long burst GRB 090926.). Parameters are listed in Tab.[3] 



larger than unity (e.g. Zhang et al. 2007). According to the standard 
internal/external shock scenario one expects the opposite, since ex- 
ternal shocks should be much more efficient than internal ones to 
dissipate the kinetic energy of the fireball. These estimates were 
based on the observed X-ray afterglow energetics (see e.g. Will- 
ingale et al. 2007; Ghisellini et al. 2009), and we can now revise 
them including the much more powerful high energy 7-ray emis- 
sion, bringing the total afterglow energetics to be roughly equal 
to the prompt phase one. Furthermore, if the fireball is indeed ra- 
diative in the first phases, with a consequent fast decay, we can 
understand why the afterglow emission at later times and at other 
frequencies is so faint. 

According to our findings, bursts detected by the LAT may be 
the ones with the largest F, and can be used to explore the high- 
end r-distribution. On the other hand, one can wonder about the 
possibility to detect with the LAT bursts with relatively smaller F, 
smaller high energy luminosities and with light curves peaking at 
larger peak times. Even if rare, nearby objects with these properties 
might be still detectable, offering a direct way to test our ideas: even 
if they should be characterised by much lower peak luminosities in 
the LAT, they should have LAT/GBM fluence ratios similar to those 
presented in this paper, and lower values of T. 

One of the argument put forward against the afterglow inter- 
pretation of the high energy flux is its variability, that according to 
Abdo et al. (2009c) can have a timescales A£ V ar as short as 90 ms. 
If true, this is certainly a severe problem for the afterglow inter- 
pretation. On the other hand the knowledge of Ai va r is limited by 
the few number of received photons. When the entire light-curve, 
lasting for a few hundreds seconds, is composed by a few hundreds 
events, one can define a very short A/ var only if there is an excep- 
tional "bunching" of photons in contiguous time-bins, and we do 
not see it in the bursts we analysed. 

Finally, we would like to emphasise the importance of estab- 
lishing, in general, if the high and low energy emission are pro- 



duced by the same electrons at the same time or instead if they are 
produced by different electrons at different times. As the study of 
GRB 0905 10 (Abdo et al. 2009b; Ghirlanda et al. 2009) has demon- 
strated, we are reaching the required data quality to put strong con- 
straints on the theories predicting the violation of the Lorentz in- 
variance at small scales, that can be tested by comparing the possi- 
ble delay of the arrival times of high energy photons. The critical 
issue about these studies is to know exactly the generation time of 
the high with respect to low energy emission. Therefore it becomes 
crucial to establish if the flux received by the LAT is the extension 
in energy of the prompt phase emission or if it is afterglow radia- 
tion. 
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